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Abstract 

The 7r-electronic structure of graphene in the presence of a modulated electric potential is in- 
vestigated by the tight-binding model. The low-energy electronic properties are strongly affected 
by the period and field strength. Such a field could modify the energy dispersions, destroy state 
degeneracy, and induce band-edge states. It should be noted that a modulated electric potential 
could make semiconducting graphene semimetallic, and that the onset period of such a transition 
relies on the field strength. There exist infinite Fermi-momentum states in sharply contrast with 
two crossing points (Dirac points) for graphene without external fields. The finite density of states 
(DOS) at the Fermi level means that there are free carriers, and, at the same time, the low DOS 
spectrum exhibits many prominent peaks, mainly owing to the band-edge states. 
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I. INTRODUCTION 



Carbon atoms could form diamond, graphites, carbon nanotubes, C60-related fullerenes, 
carbon onions, and carbon tori. These systems own special symmetric configurations, and 
their dimensionalities range from 3D to OD. These systems are of great interest to community 
because of the spectacular physical properties to which the unique geometric structures give 
rise. Recently, a new carbon-based material is discovered by the success in experimental re- 
alization of fabricating a single graphene sheet , and immediately this unique 2D system 
becomes the focus of both theoretical and the experimental researches. One of the unique 
physical properties of the kind of system, for example, is the observation of an anomalous 
quantum Hall effect [3], where the plateau of zero Hall conductance is unexpectedly absent. 
The novel phenomenon of Hall conductivity quantization is attributed to that the dynamics 
of quasiparticles in graphene is effectively relativistic, contrasting sharply with the conven- 
tional integer quantum Hall effect in the regular 2D electron system that are governed by 
the Schrodinger fermions. 

Graphene has a honeycomb crystal structure with two distinct triangular sublattices A 
and B. This real-space structure corresponds to a triangular lattice in the reciprocal space, 
which also defines a hexagonal Brioullin zone with two inequivalent corners, K and K' . The 
unique geometrical configuration of graphene leads to several non-trivial characteristics of 
its energy spectrum. There exist two identical low-energy bands in the vicinity of the two 
inequivalent corners, of which each exhibits a linear energy dispersion with conduction and 
valence bands crossing right at the corner. The highly diminished Fermi surface gives rise 
to a zero density of states at the Fermi level, so graphene a zero-gap semiconductor. 

Studying the behavior of electrons under various kinds of external fields is not only of 
vital importance for understanding their physical nature of materials but also helpful for 
designing novel devices or developing applications. Many researches has been conducted on 
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properties of graphene, such as electronic properties 
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properties, an intriguing field of less exploration [23j, |24], |2J 
graphene in the modulated fields are not yet realized experimentally; however, the 2D elec- 
tron gas formed by GaAs/AlGaAs heterojunctions in modulated fields are well-established 
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systems and have been subject of active studies in the past two decades|29|, |30|, |3l|, |32j. For 
the 2D electron gas, there are probably two types of periodic modulations in actual systems. 
One is an electrostatic modulation and the other is a magnetic-field modulation. The former 
can be realized by a periodic array of gate electrodes [30] and the latter by depositing mag- 



netic materials on the surface 
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32| . In the case of magnetic modulation, the electronic 



properties of graphene have been recently studied [27|, |28( . In this paper, we would like to 
further investigate the electronic properties of graphene in a modulated electric potential. 
The tight-binding model is employed to calculate the energy spectrum and the density of 
states. Both are analysed as functions of the field strength, the period, and the direction 
by solving the Hamiltonian matrix numerically. In this work, it is shown that semiconduct- 
ing graphene can be made semimetallic by applying a modulated electric potential. Any 
finite value of field strength will cause such a transition with a prerequisite that the period 
of modulated potential is longer enough. However, a qualitative study, recently made by 
other authors, on the Dirac particles tunneling through a one-dimensional potential bar- 
rier predicts that the 2D light-cone structure is still preserved under such a field 26[]. The 
underlying reason causing the different findings is elaborated in the discussion part. 

This paper is organized as follows. The tight-binding Hamiltonian matrix in a periodic 
electric potential is derived in Sec. II. The main characteristics of the 7r-electronic structures 
are discussed in Sec. III. Finally, Sec. IV contains concluding remarks. 



II. MODEL AND METHODS: TIGHT-BINDING MODEL 

The 7r-electronic structure of graphene is calculated by the tight-binding model within the 
nearest-neighbor interactions. Given that graphene is a two-dimensional triangular lattice 
with two-point basis, we construct the 7r-electronic eigenfunction of the system: 

l*k>=C£|^> + q*|*£), (1) 

where the tight-binding Bloch function |<&£) i s the superposition of the 2p z orbitals 

from periodic A (B) atoms. The coefficients are obtained by diagonalizing the Hamiltonian 
matrix in the space spanned by the tight-binding Bloch functions. Moreover, with the two 
prerequisites that site energies for A and B atoms are the same, and that intra-sublattice 
hoppings are ignored, we can set ( &£\Ho\&^) = ( $>^\Ho\<&^) = 0. The nonvanishing matrix 
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elements are given by 



(**\Ho\*£) 



1 



7o exp{-ik • Tj} = ^ *»k, « = 1, 2, 3, 



(2) 




where 70 = / y?p z (r — HA)H (p Pz (r — H A — Tj)d 3 r ~ 2.6 eV is the nearest-neighbor 
hopping integral with Ra being the lattice vector and ip Pz the atomic 2p z orbital. T, 
are the vectors connecting a carbon atom to its three nearest neighbors, and they are 
given by T x = (6/2, >/3&/ 2 ). T 2 = (6/2, — \/36/2), and T 3 = (-6,0), where 6 = 1.42 



given by t lk = 7 exp[ (ik x b/2 + ik y y/3b/2) ], t 2 k = 7o exp[ (ik x b/2 — ik y \/3b/2)], and 
^3k = 7o exp( — ik x b ). After specifying the matrix elements, the zero-field Hamiltonian is 
expressed as the following 2x2 Hermitian matrix, 



The Hamiltonian matrix can be solved analytically, and it generates two energy dispersions 
represented by E±(k) = ±7 [1 + 4cos 2 (V3k y b/2) + 4 cos (V3k y b/2) cos (3bk x /2)] 1 / 2 . The 
low-energy spectrum shows symmetry between two valleys (around K and K'), where both 
have the identical light-cone structure described by the relation E±(q) = ±vp\q\, with light 
velocity vf = 3670/2 and momentum q measuring the difference q = k — (k — \s.k>) in 
wavevector k and corner k^ (k^') of the Brioullin zone. 

Considering graphene exists in a ID modulated potential U(x) along the armchair di- 
rection (x), and the potential profile is assumed to take the form U(x) = Vocos(2ttx/Ie) 
with V being the field strength. For convenience, the period, l E , is further designed to 
be a multiple of 36, with l E = 3bR E - In doing so, the periodicity caused by such a field 
is made commensurate with the crystal potential of graphene itself, which helps define a 
primitive cell. The rectangular primitive cell is enlarged to include 4R E carbon atoms de- 
noted as A n (B n ) with n = 1, 2, . . . , 2R E for A(B)-type carbon atoms (Fig. 1(a)). The 
corresponding rectangular Brioullin zone shrinks to be 1/2R E of the original hexagonal 
Brioullin zone, and its dimension along the modulated direction (\k x \ < n/3bR E ) is rela- 
tively shorter than the other one (\k y \ < 7r/y / 36), as shown in Fig. 1(b). The eigenfunction 
of a system is the superpostion of elements in the basis composed of 4R E Bloch functions, 



{10k 1 ), 10k 1 ), 10k 2 ), 10k 2 ), • • • 10^"), 10?^"), l0k MB ), l0f 2flE )}, and is represented as 



A is the C-C bond length, t^, denoting individual hopping process, is then explicitly 




(3) 




(4) 
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When the period is sufficiently large, the effects of the electric potential on the off-diagonal 
matrix elements are negligible. Meanwhile, the the diagonal matrix elements would become 



($M#|$£"> =V cos[(n-l)n/R E ]=U n , 
(H>*»\H\$*>) = V cos [(n- 2/3)n/ R E ] = U n+1/3 , 



(5) 



where H = H + U is the total Hamiltonian. In effect, the assumption made above is 
plausible. The period (> 0.1 jum) in a typical experimental setup is still much longer than 
the extension of atomic 2p z orbital 1 A), so the atomic orbital feels approximately 
a constant potential within its distribution. To facilitate the calculation, the basis func- 
tions are further rearranged in a specific sequence to obtain the band Hamiltonian ma- 

trix, {|<^>, \^ 2Re ), 10k 1 ), itf"*), ■ ■ ■ \4 RE ), \C E+1 ), \C E ), \4 Re+1 )}- The band Hamil- 
tonian matrix of graphene in a modulated electric potential along the armchair direction is 
written as 



U, q p* 
Q* U 2 r e+ i/3 
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(6) 









p 

'•• U Re+1/3 q* 
p* q U Re+1 ) 

with p = ii k + ^2k and q = t 3k . 

Graphene owns the anisotropic geometry so that the 7r-electronic structure could depend 
on the modulation direction. When the modulation direction is changed into the zigzag 
direction (x), the diagonal matrix element becomes 

{^\H\^) = U n . 



k"|-f/"|$k n / 



(7) 



In this case, the period is designed as l' E = \^3bR,E for the same reason mentioned earlier in 
the armchair case. The corresponding rectangular Brioullin zone also shrinks to be 1/2Re of 
the original hexagonal Brioullin zone, and its dimension is characterized by \k x \ < n/ v^3W?e 
along the zigzag direction and \k y \ < n/3b along the other. The Hamiltonian matrix can be 
easily constructed in a similar way to that in the armchair case (not shown here). 



In the following discussion, we will consider the situation that the actual modulation 
period is about submicron length, which corresponds to Re ~ 250. This involves a process 
of diagonalizing a very large Hamiltonian matrix, and it is solved numerically to obtain the 
energy bands. Due to the resulting unoccupied conduction bands (E c ) and occupied valence 
bands (E v ) being symmetric about the Fermi level (E F = 0), we only discuss the former. 

III. ELECTRONIC PROPERTIES 

Before presenting the results under a modulated electric potential, a brief review of the 
main features of low-energy bands of the zero-field system is made. As mentioned earlier, 
the low-energy band structure has a linear energy dispersion around K or K' of the hexagon. 
However, it is not convenient to compare directly with the coming case when a modulated 
field is present, because they are presented in different Brillouin zones. To treat them on an 
equal footing, the unit cell of the zero-field system is chosen to be identical to the primitive 
cell of the system in a modulated potential. That is, if the modulated field has a period 
Re, the unit cell is chosen to be 2R E times its primitive cell. In this way, all electronic 
states in the hexagonal Brillouin zone are folded into a rectangular one. Notice that such a 
folding method guarantees not to change any features of the original band structure. Figure 
1(c) shows the fc^-dependence of the low-energy spectrum at Re = 250 for k x = (black 
solid lines) and k x = n/3bRE (blue dotted lines). There exists a nondegenerate ID linear 
band and doubly degenerate ID parabolic bands for k x = 0, while they are purely doubly 
degenerate ID parabolic bands for k x = n/3bR E - The evolution of bands from k x — to 
k x = n/3bR E fills the states between them, which reflects the 2D characteristic of the light- 
cone structure. (This energy spectrum corresponds to one of the two valleys containing 
the Dirac point K now located at k x = and k y = 2n/3^/3b.) Above the Dirac point, 
there exists a local minima for each ID parabolic band. Notice that these points should not 
be considered as band-edge states. For these states, the k x - dependent energy dispersion is 
linear with nonzero first derivative, so they cannot be treated as critical points. The band 
structure is in all respects the same as one before folding the zone. Although this alternative 
representation of band structure somewhat complicates the explanation, it is advantageous 
for the following discussion to compare the zero-field system and the system in a modulated 
electric potential. 
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A modulated electric potential has a strong effect on the energy dispersions, state degen- 
eracy, and band-edge states. The /c y -dependent conduction bands for k x = (black solid 
lines) and k x = 7r/3W?£ (blue dotted lines) are shown in Fig. 2(a) for the modulated elec- 
tric potential along the armchair direction with the period Re = 250 (~ 100 nm) and the 
strength Vq = 0.025 70. Moreover, the potential mainly affects the structure of some low ID 
bands. The original doubly degenerate parabolic bands are split. The energy dispersions 
around the k y = 2ii/3\/3b = ky are strongly deformed, and induce several band-edge states 
which are saddle points in the energy-momentum space. These band-edge states always 
appear in pairs at two sides of ky , and two band-edge states is each pair might have small 
energy difference. Far away from the ky , the spectrum is almost linear for the k y depen- 
dence, while it becomes dispersionless for the k x dependence, which can been seen in Fig. 
2(a) where the bands for k x = and k x = 7r/3W?£ have identical dispersion. It is important 
to note that there exist more Fermi momentum states fcp's. The energy dispersions near 
fc^'s are linear for the /c^-dependence, while they are completely flat for the /^-dependence. 
The dispersionless feature means that the number of the Fermi-momentum states is infinite, 
which sharply contrasts with the original two Fermi-momentum states or Dirac points in 
zero-field graphene. Besides, the non-zero measure of the Fermi surface indicates finite value 
of the DOS at the Fermi level. 

The aforementioned features might rely on the strength, period, and direction of a modu- 
lated electric potential. The modulation effects are even enhanced when field is strengthened, 
illustrated in the Fig. 2(b). When the strength increases, more low ID bands are affected. 
The energy bands are further deformed, and more band-edge states are created. In addition, 
the number of Fermi-momentum states increases almost linearly with the field strength. For 
the case in Fig. 2(b), the number of Fermi-momentum states is approximately double that 
in Fig. 2(a). The modulation period can also affect the numbers of band-edge states and 
Fermi-momentum states, as shown in Fig. 2(c). When the period increases, more band-edge 
states are created. The number of Fermi- momentum states remains almost unchanged in 
the long-period regime, whereas there is a drastic change in the small-period regime. The 
former can be understood by taking into consideration both the accompanying reduced k x 
length and increased fop's associated with different kyS in accordance with the increasing 
period, while the latter will be discussed later by calculating the density of states at the 
Fermi level. The low-energy bands are, on the other hand, not sensitive to the direction of a 
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modulated electric potential. Figure 2(d), for example, presents the bands for the potential 
modulated along the zigzag direction. 

The essential features of electronic structure are directly reflected on the density of states, 
and defined as 



Because of the linear energy dispersion, the low-frequency DOS without fields is proportional 
to iv, as shown by dotted line in Fig. 3(a), and it has no special structures. The vanishing 
DOS at Ep = indicates that graphene is a zero-gap semiconductor. A modulated electric 
potential can lead to many prominent peaks and a finite DOS at uj = (Figs. 3(a)-3(d)). 
The peak structures come from the band-edge states of parabolic bands (Figs. 2(a)-2(d)). 
The frequency, number, and height of peaks are sensitive to the changes in the field strength 
and period. The DOS at the Fermi level means that there are free carriers, so graphene 
becomes a semimetal in the presence of a modulated electric field. The value of DOS at 
iv = grows as the field strength increases (Figs. 3(a) and 3(b)). However, it does not 
change significantly as the period is varied (Fig. 3(c)). Moreover, the low-frequency DOS 
is unaltered when the modulation direction changes (Fig. 3(d)), which reflects the isotropic 
symmetry of graphene in low-energy spectrum. 

The singularities in the DOS might cause special structures in some physical quantities, 
for example, giving rise to the strong absorption peaks in optical measurement, so it is 
worthy to investigate their properties in detail. The peak height, number and frequency are 
dominated by the strength and period of the potential (Figs. 3(a)-3(c)). The peak height 
is enhanced by the increasing strength, whereas it is reduced by the increasing period. The 
peak number is increased by both the increasing strength and period. The relations between 
of the peak frequency (ub e ) and the field condition are elaborated through examining the 
first four prominent peaks, as shown in the Fig. 4(a) for the strength dependence and Fig. 
4(b) for the period dependence. From these, it is observed that the peak frequency weakly 
depends on the strength, while declines and presents somewhat oscillatory behavior as the 
period increases. 

The relationship between the DOS at uj = and the field strength and that between the 
field strength and the period deserve a closer investigation. Referring to Fig. 5(a), D(uo = 0) 
shows a nearly linear variation with the field strength. This relation just reflects the fact 
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that the number of Fermi-momentum states is approximately proportional to field strength. 
On the other hand, the period dependence of D(u = 0) exhibits different features in short- 
and long-period regimes, as shown in Fig. 5(b). The two regions are distinguished by a 
threshold period R t h, which happens at Re ~ 100 for Vq = 0.025 70 ( black triangles in Fig. 
5(b)) and becomes shorter as the field strength gets stronger ( red circles in Fig. 5(b)). For 
Re -C Rth, D(u = 0) remains zero as Re increases from zero, which implies that the 2D 
light-cone structure is not affected at such a short period. As Re is increased to cross R t h, 
D(uj = 0) quickly grows and saturates to a definite value, with some fluctuations in the long- 
period regime. If one consider a realistic experimental system, the typical modulation period 
will be about submicron or even longer. The D(u = 0) then would fall into the long-period 
regime, so its value does not change significantly as the period is varied. Therefore, this 
feature should be very robust even there exist some inevitable period fluctuations coming 
from processing devices or experimental setups. 

These interesting features, including field induced band-edge states and the semimetallic 
behavior of graphene, caused by the modulated electric potential are the main results in this 
work. Nevertheless, the semimetallic behavior, the most important feature, was not obtained 
in a recent study The study considered graphene subjected to a square- wave like 

modulated potential, concluding that the 2D light-cone structure is still preserved. Although 
the conclusion seems to contradict with our results given here, a in-depth examination shows 
that this is not the case. The main discrepancy lies in the fact that the modulation period 
used (20 nm pa Re = 50) in their calculation is too short to achieve the threshold {Rth), 
so that the semimetallic behavior was not observed. To clarify this point, the dependence 
of D(u = 0) on the field strength and the modulation period for a square- wave potential is 
calculated here, as shown by green squares in Figs. 5(a) and 5(b) respectively. D(uj = 0) 
in Fig. 5(a) exhibits a linear relation with the field strength, which is similar to the cosine- 
wave potential. The threshold period in this case is Rth ~ 100 (Fig. 5(b)), which is close to 
that in the cosine- wave potential. As a result, from Fig. 5(b), it self-evident ly explains the 
reason why the semimetallic behavior is not observed by other authors solely because that 
Re taken in their calculation is less than R t h- 
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IV. CONCLUDING REMARKS 



In summary, the tight-binding model is used to investigate the effects of a modulated 
electric potential on the 7r-electronic structure of graphene. The low-energy electronic prop- 
erties are mainly dominated by field strength and period. However, they are not sensitive 
to the field direction, which is due to the isotropy of the original low-energy bands. A mod- 
ulated electric potential drastically changes the energy dispersions, state degeneracy, and 
band-edge states. The density of states exhibits a lot of prominent peaks. These structures 
are mainly determined by the energy dispersions and the band-edge states. Notice that 
there are infinite Fermi-momentum states. The finite DOS at the Fermi level indicates the 
existence of free carriers. The semiconducting graphene becomes semimetallic by applying a 
modulated electric potential, as long as the period surpasses a threshold period (Re 3> Rth)- 
Different nature of the low-energy excitations for graphene in the presence of a modulated 
potential is expected to give rise intriguing physical properties. For instance, the free elec- 
trons are expected to cause the low-frequency plasmon. The theoretical predictions could 
be tested by the experimental measurements on the energy loss spectra. 
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FIG. 1: (a) The rectangular primitive cell of graphene in a modulated electric potential along the 
armchair direction (x). The rectangular first Brillouin zone is shown in (b), and the hexagon is 
that without external fields. The /c^-dependent low-energy bands without any fields at Re = 250 
is shown in (c) for k x = (black solid lines) and k x = 27r/3bRE (blue dotted lines). 



13 



0.08 



0.06 - 



o 



0.04 - 



0.02 - 



\w\ww\'- \ — 7 TTimum 

(a) armchair, R ) =250; V =0.025 7„ 





(b) armchair, R =250; V =0.05 7, 




0.65 0.7 

fa ( n /^b) 



0.65 0.7 
fcy (VV3b) 



/.■■■//■Ill 




0.65 0.7 -0.05 

ky (VV3&) ^ (V3fc) 



0.05 



FIG. 2: The /c^-dependent low-energy bands in the modulated field along armchair direction with 
Vq = 0.025 70 and Re = 250 for k x = (black solid lines) and k x = 2ir/3bRE (blue dotted lines), 
those in a stronger field with Vq = 0.05 70, and those in the field of a longer period Re = 500 are 
respectively shown in (a), (b), and (c). (d) is the ^-dependence of low-energy bands at Vq = 0.025 
70, Re = 433 with the potential modulating along the zigzag direction for k x = (black solid 
lines) and k x = 2tt / y/3bRE (blue dotted lines). Notice that the actual periods in (a) and (d) are 
almost of equal length, and that the units of k y are different for armchair (-7r/v / 3ft) and zigzag 
[it /3b) direction. 14 
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FIG. 3: The low- frequency density of states for the potential modulated along armchair direction 
at (a) V = 0.025 70, R E = 250; (b) V = 0.05 70, R E = 250; (c) V = 0.025 70, R E = 500; (d) 
along the zigzag direction at (Vb = 0.025 70, Re = 433). The dashed line in (a) is for the zero-field 
spectrum. 
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FIG. 4: Energies of the first four prominent peaks versus (a) the field strength at Re = 500 and 
(b) the modulation period at Vq = 0.05 70. 
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FIG. 5: The density of states at the Fermi level versus (a) the field strength at Re = 250 and 
500 (Re = 250) for the modulated potential of cosine- wave (square- wave) profile, and versus (b) 
the modulation period at Vo = 0.025 and 0.05 70 (Vq = 0.025 70) for that of the cosine-wave 
(square- wave) profile. 
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